Confronting Dark Energy Models using Galaxy Cluster Number Counts 

S. BasilakofQ 

Academy of Athens, Research Center for Astronomy and Applied Mathematics, Soranou Efesiou 4, 11527, Athens, Greece 

M. PlionisQ 

Institute of Astronomy & Astrophysics, Nationals Observatory of Athens, Thessio 11810, Athens, Greece, and 
Instituto Nacional de Astrofisica, Optica y Electronica, 72000 Puebla, Mexico 

J. A. S. LimM 

Departamento de Astronomia (IAGUSP), Universidade de Sao Paulo 
Rua do Matao, 1226, 05508-900, S. Paulo, Brazil 

The mass function of cluster-size halos and their redshift distribution are computed for 12 distinct 
accelerating cosmological scenarios and confronted to the predictions of the conventional flat ACDM 
model. The comparison with ACDM is performed by a two-step process. Firstly, we determine the 
free parameters of all models through a joint analysis involving the latest cosmological data, using 
SNe type la, the CMB shift parameter and BAO. Apart from a brane world inspired cosmology, 
it is found that the derived Hubble relation of the remaining models reproduce the ACDM results 
approximately with the same degree of statistical confidence. Secondly, in order to attempt distin- 
guish the different dark energy models from the expectations of ACDM, we analyze the predicted 
cluster-size halo redshift distribution on the basis of two future cluster surveys: (i) an X-ray sur- 
vey based on the eROSITA satellite, and (ii) a Sunayev-Zeldovich survey based on the South Pole 
Telescope. As a result, we find that the predictions of 8 out of 12 dark energy models can be 
clearly distinguished from the ACDM cosmology, while the predictions of 4 models are statistically 
equivalent to those of the ACDM model, as far as the expected cluster mass function and redshift 
distribution are concerned. The present analysis suggest that such a technique appears to be very 
competitive to independent tests probing the late time evolution of the Universe and the associated 
dark energy effects. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36.-fx 



1. INTRODUCTION 

Recent studies in observational cosmology lead to the 
conclusion that the available high quality cosmological 
data, from Supernovae type la (SNe la), matter power 
spectrum analysis, the angular power spectrum of the 
cosmic microwave background (CMB), baryon acoustic 
oscillations (BAO) and other complementary probes, like 
the existence of old galaxies at high redshifts, are well fit- 
ted by an emerging "standard cosmological model" (see 
[IHH and references therein). By assuming flatness, as 
theoretically required by inflation and the CMB observa- 
tions, the standard model is described by the Friedmann 
equation: 

H 2 (a) = Q = ^ \p m {a) + pdb(o)] , (1.1) 

where a(t) is the scale factor of the universe, p m (a) is the 
density corresponding to the sum of baryonic plus cold 
dark matter (the latter needed to explain clustering) , and 
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an extra component /0de(oO; with negative pressure called 
dark energy (DE hereafter), that accelerates the cosmic 
expansion (for reviews see @). 

Nowadays, the physics of dark energy is considered one 
of the most fundamental and challenging problems on 
the interface uniting Astronomy, Cosmology and Parti- 
cle Physics. In the last decade there have been theoret- 
ical debates among the cosmologists regarding the na- 
ture of this exotic component. Many candidates have 
been proposed in the literature, such as a cosmolog- 
ical constant A (vacuum), time- varying A(t) cosmolo- 
gies, quintessence, k— essence, quartessence, vector fields, 
phantom, tachyons, modifications ofgravity, Chaply- 
gin gas and the list goes on (see [lfj| - r26j and references 
therein). Generically, some high energy field theories also 
suggest that the equation of state of such a dark energy 
is a function of cosmic time (see, for instance, .27] ). Nat- 
urally, in order to establish the evolution of its equation 
of state (EoS), a realistic form of H(a) is required and 
must be constrained through a combination of indepen- 
dent DE probes. 

On the other hand, the abundance of collapsed struc- 
tures as a function of mass and redshift is a key statisti- 
cal test for studies of matter distribution in the universe, 
and, more importantly, it can readily be accessed from 
observations [28(. Indeed, the mass function of galaxy 
clusters has been measured based on X-ray surveys [291 — 
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|31|, via weak and strong lensing studies |32| - |34| |, using op- 
tical surveys, like the SDSS I3& 1391 . as well as, through 
Sunayev-Zeldovich (SZ) effect [37|. In the last decade 
many authors have been involved in this kind of stud- 
ies and have found that the abundance of the collapsed 
structures is affected by the presence of a dark energy 
component [38l - l46j . 

In this work, we discuss how to differentiate a large 
family of flat DE cosmologies (12 models) from the con- 
ventional ACDM model. Initially, a joint statistical anal- 
ysis, involving the latest observational data from SNe 
type la, CMB shift parameter and BAO, is implemented. 
Since however the resulting Hubble functions of these DE 
models (apart a brane world scenario) are statistically in- 
distinguishable, we attempt to discriminate the different 
DE models by computing the halo mass function and the 
corresponding redshift distribution of the cluster-size ha- 
los. Finally, by using future X-ray and SZ surveys we 
show that the evolution of the cluster abundances (espe- 
cially at large redshifts 0.5), is a potential discrimina- 
tor (in the Popper sense) for a large fraction of the stud- 
ied DE models. As an extra bonus, we find that many 
of the DE models could be differentiated, from ACDM, 
even with the present available observational data from 
galaxy cluster mass function. 

The article is planned as follows. Section 2 is a pro- 
vision of dark energy issue and the overall approach 
adopted in the paper. In section 3, a joint statistical 
analysis based on SNe la, CMB and BAO is carried out 
for the conventional ACDM model. This statistical ap- 
proach is used to constrain the DE model parameters and 
it is presented in section 4. The linear growth factor of all 
models is discussed in section 5, while in section 6, we dis- 
cuss and compare the corresponding theoretical predic- 
tions regarding the evolution of the cluster abundances. 
Finally, the main conclusions are summarized in section 
7. Throughout the paper we adopt Hq = 71km/sec/Mpc. 



EoS parameter in terms of E(a) — H{cl)/Hq [13] using 
the Friedmann equations as 
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Since the exact nature of the DE has yet to be found, 
the above DE EoS parameter encodes our ignorance re- 
garding the physical mechanism powering the late time 
cosmic acceleration. 

The methodology described above can also be applied 
to the framework of modified gravity (see [23|, l48j|). In 
this case, instead of using the exact Hubble flow through 
a modification of the Friedmann equation one may con- 
sider an equivalent Hubble flow somewhat mimicking Eq. 
(|1.1[) . The key point here is that the accelerating ex- 
pansion can be attributed to a kind of "geometrical" 
DE contribution. Now, since the matter density (bary- 
onic+dark) cannot accelerate the cosmic expansion, we 
perform the following parametrization [23l |48| : 



E 2 (a) 



H 2 (a) 
Hi 



5H 2 



(2.3) 



Naturally, any modification to the Friedmann equation 
of general relativity may be included in the last term of 
the above expression. After some algebra one may also 
derive, using Eqs. (|2.2p and (|2.3p . an effective ("geomet- 
rical") dark energy EoS parameter, given by: 



w(a)- 1 1 dlnSH2 
3 dlna 



(2.4) 



The above formulation will be adopted in our statistical 
analysis of all DE models discussed in section 4. 



3. LIKELIHOOD ANALYSIS: THE ACDM CASE 



2. DARK ENERGY EQUATION OF STATE 

In what follows, it will be assumed that the universe 
is a self-gravitating fluid described by general relativ- 
ity, and endowed with a spatially flat homogeneous and 
isotropic geometry. In addition, we also consider that 
it is filled by non-relativistic matter plus a DE compo- 
nent (or some effective mechanism that simulates it), and 
whose equation of state, pde = w(a)/9DE, is driving the 
present accelerating stage. Following standard lines, the 
Hubble parameter reads: 



E 2 (a) 



H 2 (a) 
Hi 
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3 fUlny[l+w(y)] 



(2.1) 



where E(a) is the normalized Hubble flow, Q m is the di- 
mensionlcss matter density at the present epoch, Qde = 
1 — f2 m denotes the DE density parameter, and w(a) its 
EoS parameter. On the other hand, we can express the 



In this section we briefly present the basic observa- 
tional samples and statistical analysis tools that will be 
used to constrain the cosmological parameters of the DE 
models. Here we discuss the ACDM model since it is now 
widely believed that if a better cosmology is called for, 
it will describe a cosmos that looks much like a ACDM 
model. 

(i) Supernovas. In our statistical analysis we consider 
the Constitution set, containing 397 SNe type la, as com- 
piled by Hicken et al. 6]. In order to avoid possible 
problems related to the local bulk flow, we use a subset 
of this sample containing 366 SNe la all with redshifts 
z > 0.02. The likelihood estimator is determined by a 
XsNia statistics 
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XsNIa(P) = 



where aj = (1 
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(3.1) 



is the scale factor of the Universe 



in the observed redshift Zi, /i is the distance modulus 
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fi = m — M = 51ogffi + 25 and d^ia, p) is the luminosity 
distance, d L {a,p) = ca- 1 J a y ^ J {y ) ■ 

(ii) Baryon Acoustic Oscillations (BAO). In addition 
to the SNe la data, we also consider the BAO scale 
produced in the last scattering surface by the competi- 
tion between the pressure of the coupled baryon-photon 
fluid and gravity. The resulting acoustic waves leaves 
(in the course of the evolution) an overdensity signature 
at certain length scales of the matter distribution. Ev- 
idence of this excess was recently found in the cluster- 
ing properties of the SDSS galaxies [49| . and it provides 
a suitable "standard ruler" for constraining DE mod- 
els. In particular, we consider the following estimator 

-,2/3 



SDSS data to be A 
a s = (1 + Zs)- 1 ~ 



I 1 ~ 

J a s a 



measured from the 
0.35 [or 



0.469 ± 0.017, where z s 

Therefore, the corresponding 



0.75] 



Xbao 



function can be written as 



Xbao(p) = 



L4(p) - 0.469] 5 
0.017 2 



(3.2) 



where p is a vector containing the cosmological fitting 
parameters. 

(iii) CMB Shift Parameter. Another interesting geo- 
metrical probe for dark energy is provided by the an- 
gular scale of the sound horizon at the last scattering 
surface. It is encoded in the location lf T of the first 
peak of the angular (CMB) power spectrum [5(| [HJ, and 
may be defined by the quantity 7Z = \JVl m a -± d E(a) ■ 
The shift parameter measured from the WMAP 7-years 
data [| is K = 1.726 ± 0.019 at z ls = 1091.36 [or 
ai s = (1 + zis)^ 1 ~ 9.154 x 10 -4 ]. In this case, the 

Xcmb 



xLb(p) 



[R(p) - 1.726] 2 



(3.3) 



0.019 2 

It should be stressed that for CMB shift parameter, the 
contribution of the radiative component, (f2^a -4 , where 
Qr ~ 4.174 x 10 -5 ft. -2 ) needs also to be considered [5j. 
Note also that the measured CMB shift parameter is 
somewhat model dependent but mostly to models which 
are not included in our analysis. For example, such is the 
case when massive neutrinos are included or when there 
is a strongly time varying equation of state parameter. 
The robustness of the shift parameter has been tested 
and discussed in 

In order to put tighter constraints on the correspond- 
ing parameter space of each DE model, the probes de- 
scribed above must be combined through a joint like- 
lihood analysis 1 , given by the product of the individual 



likelihoods according to: C to t{p) = £sNia x £bao x Ccmb, 
which translates in the joint x 2 function in an addition: 
XL(P) = XsNia + Xbao + Xcmb- The results of this anal- 
ysis provide: Q m = 0.28 ± 0.01 with xf ot (O m )/<io/ ~ 
439.6/367. 

In this concern, it should be remarked that a value of 
XsNia — 439 is what one gets by using directly the Con- 
stitution SNIa set of Hicken et al. [6J since the number of 
the SNIa used (366) dominates the overall Xtot budget of 
the joint likelihood. In addition, the corresponding good- 
ness of fit {x 2 SNia/d°f — 1-2) is significantly larger then 
the one found by Davis et al. Q {XsNia/dof — 1) from a 
previous sample containing 192 SNIa (see the discussion 
in [HI, HH ) . Such a discrepancy appears to be the out- 
come of the different approaches chosen in order to join 
the different contributing SNIa sets. According to Hicken 
(private communication, 2009) in the case of the Davis 
et al. 0] data the nearby SNIa were imposed to provide 
X 2 1 1 dof ~ 1 by hand, while no such fine-tuning was im- 
posed on the Constitution set. Also the latter set includes 
distant SNIa which have typically larger distance modu- 
lus uncertainties with respect to those used in Davis et 
al. sample 0- In particular, this means that the higher 
XsNia/dof value of the Constitution set should probably 
be attributed to typically lower SNIa distance modulus 
uncertainties. 2 

By using the most recent BAO results of Percival et 
al. [55[, we have checked and verified that the above 
constraints of the flat ACDM model are not apprecia- 
bly modified. Therefore, the reader should keep in mind 
the value Xtot — 439.6 for any further comparison with 
results predicted by the alternative DE models. 



4. CONSTRAINTS ON DARK ENERGY 
MODELS 

Let us now present the twelve flat DE models whose 
free parameters will be constrained by using the same 
methodology and cosmological data as that applied to 
the ACDM model (see previous section). 



A. Constant Equation of State (XCDM model) 

In this kind of cosmology (hereafter XCDM-models) 
the equation of state parameter is constant [56]. In a 
point of fact, these DE models do not have a clear phys- 
ical motivation. In particular, for quintessence models 
driven by a real scalar field, a constant EoS parameter 
requires an extreme fine tuning of its potential and ki- 
netic energy. In spite of that, this subclass of DE models 



Likelihoods are normalized to their maximum values. In the 
present analysis we always report la uncertainties on the fitted 
parameters. Note also that the total number of data points used 
here is Ntot = 368, while the associated degrees of freedom are: 
dof= Ntot — iflt, with nut the model-dependent number of fitted 
parameters. 



2 Such a possibility has been crudely tested by Plionis et al. J54j. A 
lower value X%_\ i a /d°f — 1-07 (similar to that of Davies 01) has 
been obtained by increasing the distance modulus uncertainty of 
the Constitution nearby SNIa (z < 0.4) by 20%. 
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have been widely used in the literature due to their sim- 
plicity. Notice that DE models with a canonical kinetic 
term have — 1 < w < —1/3. Models with (w < — 1), 
sometimes called phantom dark energy [57j | , are endowed 
with a very exotic nature, like a scalar field with negative 
kinetic energy. Now, by using Eq. (|2.1j) the normalized 
Hubble parameter becomes 



E 2 {a) = ft r , 



d-n r , 



-3(l+tu) 



(4.1) 



Now, in order to constrain XCDM-models with the 
observational data we sample VL m € [0.1,1] and w € 
[—2, —0.4] in steps of 0.01. As a result, we find that the 
best fit values are Q m = 0.28±0.01 and w = -0.99±0.05 
with Xtot(tt m ,w)/dof ~ 439.5/366. 

The above results are in excellent agreement with those 
found by different authors 0, 0-0, HI, 31 1 . It is also worth 
noticing that the concordance ACDM cosmology is de- 
scribed by a XCDM model with w strictly equal to -1. 
The corresponding limits form the basis of our present 
comparison and were separately derived (see last section) . 



B. Parametric Dark Energy (CPL model) 

This kind of cosmology was first discussed with basis 
on the Chevalier-Polarski-Lindcr [58, 59] parametrization 
(hereafter CPL). The dark energy EoS parameter is de- 
fined as a first order Taylor expansion around the present 
epoch: 



w(a) = wq + wi(l — a), 



(4.2) 



where wq and w\ are constants. The normalized Hubble 
parameter now reads: 



E 2 (a) = O r 



-(1-0, 



-3(l+«)o+i"i) e 3Mii(a-l) 



(4.3) 



functional form of the normalized Hubble flow as given 
by Eq. (|2.3|) is affected only by one free parameter, fi m . 
The quantity SH 2 is given by 



SH 2 = 2il h 



(4.4) 



where Qb w — (1 — O m ) 2 /4. From Eq |2.4[ it is readily 
checked that the geometrical DE equation of state pa- 
rameter reduces to 



w(a) = 



1 



1 + £l m (a) 



(4.5) 



where f2 m (a) = £l m a~ 3 / E 2 (a) . The joint likelihood anal- 
ysis provides a best fit value of fl m = 0.23 ±0.01, but the 
fit is much worse, xiot(^m)/dof ~ 500/367, in compari- 
son with the one provided by a ACDM cosmology. 



D. Low Ricci Dark Energy (LRDE model) 

In this model we use a simple parametrization for the 
Ricci scalar which is based on a Taylor expansion around 
the present time: TZ(a) = tq +t*i(1 — a). It is interesting 
to mention that at the early epochs the cosmic evolution 
tends asymptotically to be matter dominated (for more 
details see [231] )■ In this framework, we have 



E 2 (a) 



a 4(r +ri-l) e 4ri(l-a) 
O m (J 



a > a t 
a < a t 



(4.6) 



where at — 1 — (1 — 4ro)/4ri. At present, the matter 
density parameter is directly related with the above con- 
stants via 



a. 



4r + 4ri 
4ri 



1 



4ro+4ri — 1 



3 l-4r 



The corresponding EoS parameter is given by 



In order to constrain the free parameters, we sample 
them as follows: w G [-2, -0.4] and w\ G [-2.6, 2.6] w ( a ) 
i n steps of 0.01. By fixing a prior in the density parame- 
ter, fi m = 0.28, we find that the overall likelihood func- 
tion peaks at w = -0.96 ± 0.13 and wi = -0.40 ± 0.70 
in very good agreement with the 7 years WMAP data 
The corresponding Xtot( w o,wi)/dof is 439.5/366. A 
value that should be compared with the ACDM predic- 
tion. 
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(4.7) 



(4.8) 



Notice, that we sample the unknown parameters as fol- 
lows: r € [0.5,1.5] and r\ € [—2.4,-0.1] and here are 
the results: r = 0.80 ± 0.02 and r x = -0.69 ± 0.05 
(fi m ~ 0.27) with Xtot(ro,n)/dof ~ 439.9/366. 



E. High Ricci Dark Energy (HRDE model) 



C. Braneworld Gravity (BRG model) 

In the context of a braneworld cosmology (hereafter 
BRG) the accelerated expansion of the universe can be 
explained by a modification of the gravitational interac- 
tion in which gravity itself becomes weak at very large 
distances (close to the Hubble scale) due to the fact that 
our four dimensional brane survives into an extra dimen- 
sional manifold (see [60| and references therein). The in- 
teresting aspect of this scenario is that the corresponding 



A different Ricci DE model was proposed by Linder 
61]. In this framework, the Ricci scalar evolves at high 
redshifts obeying the following expression 



1 — 3wq- 



SH 2 



H 2 



(4.9) 



where SH 2 = E 2 (a) — H m a 3 . In this model the normal- 
ized Hubble flow becomes: 



E 2 (a) = n m a~ 3 (l + Pa^ W0 ) 



3uM -lnQ m /ln(l+/3) 



(4.10) 
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where f3 — O^ 1 — 1. As in the previous case, the effec- 
tive equation of state parameter is again related to E(a) 
according to Eq. (|2.2[) . Now, by using the same sampling 
as in the XCDM-models, we find that the joint likelihood 
peaks at O m = 0.28 ± 0.01 and w = -0.99 ± 0.05 with 
x 2 tot (n m ,w )/dof ~ 439.5/366. 



F. Tension of Cosmic Magnetic Fields 
(TCM model) 

A couple of years ago, Contoupolos and Basilakos [62| 
proposed a novel idea which is based on the following con- 
sideration (hereafter TCM): if the cosmic magnetic field 
is generated in sources (such as galaxy clusters) whose 
overall dimensions remain unchanged during the expan- 
sion of the Universe, the stretching of this field by the 
cosmic expansion generates a tension (negative pressure) 
that could possibly explain a small fraction of the DE 
(~ 2 — 5%). In this flat cosmological scenario the nor- 
malized Hubble flow becomes: 



E 2 {a) = n m a' 3 + SH 2 , SH 2 = A + O s a 



-3+n 



(4.11) 



where £Ib is the density parameter for the cosmic mag- 
netic fields and fl\ = 1 — O m — From Eq. (|2.4I) it 
is readily seen that the effective EoS parameter which is 
related to magnetic tension reads 



w(a) 



30a + nflga 
3(Oa + O^a 



-3+n 
-3+n\ 



(4.12) 



Again, in order to constrain the parameters, we sample 
Qb € [0.,0.3] and n € [0, 10]. By considering a prior of 
TO = 0.28, the best fit values are: O s = 0.01 ± 0.005 
and n = 1.80 ± 0.80 with X 2 t(^B,n)/dof ~ 439.8/366. 



G. Pseudo-Nambu Goldstone Boson 
(PNGB model) 

In this cosmological model the DE equation of state 
parameter is expressed with the aid of the potential 
V{cj)) cx [l + cos(<f>/F)] Hi: 



w(a) = -1 + (I+Wq)^ 



(4.13) 



where w <E [-2,-0.4] and F £ [0,8]. We see that for 
n«l, the EoS parameter goes to w(a) = —1. Based on 
this parametrization the quantity SH 2 (see eq |2.4[) takes 
the following form 



SH 2 = (1 - O m )exp 



3(1+ wp) 
F 



(1 - a F ) 



(4.14) 



In this case, our joint statistical analysis yields that the 
likelihood function peaks at w = —0.97+ 0.09 and F — 
5.9 ± 3.2 with x 2 ot( w o, F)/dof ~ 439.7/366. 



H. Early Dark Energy (EDE model) 

Another cosmological scenario that we include in our 
paper is the early dark energy model (hereafter EDE). 
The basic assumption here is that at early epochs the 
amount of DE is not negligible 64|. In this framework, 
the overall DE component is given by 



Ode (a) 



1-fir. 



n e (i 



-3wo 



n e (i~ a - 3wo ) 



i - n m + n m a 3w ° 

where O e is the early DE density and wq is the equation 
of state parameter at the present epoch. We observe that 
the EDE model was designed to satisfy the two following 
conditions: (i) mimic the effects of a late dark energy, 
and (ii) provide a decelerating expansion of the Universe 
at high redshifts. The normalized Hubble parameter is 
now written as 



E 2 (a) = 



-3 



1 - DE (a) 



(4.16) 



and from Eq. (|2.2p . one may obtain the EoS parameter 
as a function of the scale factor. 

Now, from the joint likelihood analysis we find that 



O e = 0.02 ± 0.02 and w = -1.04 



+0.07 
-0.11 



(for the prior 



of Q m = 0.28) with xiot(^e,w )/dof ~ 439.3/366. Note, 
that a different class of EDE models was recently studied 
in El. 



I. Variable Chaplygin Gas (VCG model) 

Accelerating cosmologies can also be driven by the so 
called variable Chaplygin gas (hereafter VCG) which cor- 
responds to a Born-Infeld tachyon action (2J, [65| . Re- 
cently, an interesting family of Chaplygin gas models 
was found to be consistent with the current observational 
data [66[ . In the framework of a spatially flat geometry, it 
can be shown that the normalized Hubble function takes 
the following formula: 



E 2 {a) = n b a- 3 + (l-n b )^B s a-s + {1 - B s 



(4.17) 



where Of, ~ 0.0226/i~ 2 is the density parameter for the 
baryonic matter [j| and B s G [0.01, 0.51] and n e [-4, 4]. 
The effective matter density parameter is: O m = Of, + 
(1 — ilb)\/B~ s - We find that the best fit parameters are 
B s = 0.05 ± 0.01 and n = 1.65^'^ (fi m ~ 0.26) with 
X 2 tot (B s> n)/dof ~ 441.3/366. 



J. Time Varying Vacuum (Arc and Aps models) 

Let us now consider the possibility of a decaying A- 
cosmology, that is, A = A(a). In this kind of model, the 
value of A is small because the Universe is too old, and, 
therefore, it alleviates the cosmological constant prob- 
lem. As it appears, the decaying vacuum equation of 
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state does not depend on whether A is strictly constant 
or variable. Therefore, the EoS takes the usual form, 
P A (t) = -p A (t) = -A(*)/8ttG. The global dynamics of 
such models have been investigated extensively in the lit- 
erature, in fact, much before the discovery of the present 
accelerating stage [l3| (for a recent and extensive dis- 
cussion see Basilakos et al. 14]). By introducing in the 
Friedman equations the idea of a time-dependent vacuum 
one obtains 

// + -//' = -. (4.18) 



The traditional A = const cosmology can be described 
directly by the integration of Eq. (|4. 1 8[) . but the same 
equation is also valid for A = A(a). This means that 
a supplementary equation for the time evolution of A is 
needed in order to unveil the dynamics of this pattern. In 
this work we consider two different versions of the A (a) 
models, namely renormalization group of the quantum 
field vacuum (2f| [63] (hereafter Arg) and a power series 
vacuum [68| (hereafter Aps) 



A(H) 



A + 3 7 (H 2 - H 2 ) for A RG 
H {j-3n m )H + (3- 7 )ff 2 forAp S 

Naturally, the vacuum energy density is normalized to its 
present value: A = A(H ) = 3Q\Hq (Qa = 1 — ^m)- 
Inserting the above expressions into Eq. (|4.18p we finally 
obtain the normalized Hubble flow 



E(a) 



1—Q„ 



1-7 
1 - 



3Q„ 
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1-7 



,-3(1-7) 



3n„ 

7 



1/2 



,-7/2 



for Arg 
for Aps 

(4.20) 

It has been shown [67|,[20] that for A(t) there is a cou- 
pling between the time-dependent vacuum and matter 
component. Indeed, by combining the conservation of 
the total energy with the variation of the vacuum en- 
ergy, one can shown that the A(t) models provide either 
a particle production process or that the mass of the dark 
matter particles increases [69| . 

Now, by applying our statistical procedure for both 
cosmologies, the best fit parameters are: 

Arg model: il m = 0.28 ± 0.01 and 7 = 0.002 ± 0.001 
with Xtot(n m ,7)/dof ~ 439.5/366. 

A PS model: ft m = 0.32 ± 0.02 and 7 = 3.45±{^ with 
X 2 ot (0 m)7 )/d O /~ 440.7/366. 



K. Creation Cold Dark Matter (CCDM Model) 



dark matter particles), and the energy conservation laws 
for each component have been investigated thoroughly by 
[72l - l74j ] . In this framework, the basic equation governing 
the global dynamics of a Universe endowed with a flat 
geometry is |72T475l | 



H 



-H z = 



3 H' 



(4.21) 



where pdm is the dark matter energy densities and T is 
the so called creation rate of the cold dark matter and 
it has units of (time) -1 . The creation pressure is nega- 
tive and defined in terms of the creation rate and other 
cosmological quantities by the expression (72Tj75| 



Pc 



PdrnF 
3H 



(4.22) 



In CCDM models, the functional form of T is phenomeno- 
logically parametrized [74j by the following relation 



r = zn A 



Pc 



(4.23) 



where f^A (called a in the 74]) is a constant, p co — 
3Hq/8wG is the present day value of the critical den- 
sity. Performing now the integration of Eq. (|4.21l) we 
obtain 



E 2 (a) = Q, 



n A 



(4.24) 



where Cl A — 1 — Ct m . Interestingly, the CCDM model re- 
sembles the global dynamics of the concordance ACDM 
cosmology by including only one free parameter. In 
particular, this means that Cl m — 0.28 ± 0.01 with 
x 2 tot (n m )/dof ~ 439.6/367. 



5. EVOLUTION OF MATTER 
PERTURBATIONS 

In this section we generalize the basic equations which 
govern the behavior of the matter perturbations within 
the framework of the previously described DE models. 
We focus our attention to the most generic models, ie., 
the A(t) and CCDM models. In these cases, a neo- 
Newtonian description can be introduced by considering 
an extended continuity equation together with the Euler 
and Poisson equations (7g, l79j |. In virtue of the parti- 
cle creation model, the evolution equation of the growth 
factor becomes [lj, [75| : 



In the context of the standard general relativity the- 
ory, is also possible to reduce the dark sector by consid- 
ering the presence of the gravitationally induced parti- 
cle creation mechanism. Like A(t)-models this kind of 
scenario has also been discussed long ago [7lj]. More re- 
cently, the cosmological equations for the mixture of ra- 
diation, baryons and cold dark matter (with creation of 



D + (2H + Q)t) - UnGpm - 2HQ - Qj D = 0, (5.1) 

where p m is the matter density, G = G and 

0(f \ _ f -A/p m for A RG or A PS , , 

\ 3Q APco H/p m for CCDM V ' ' 



7 



In the case of non interacting DE models, [Q(t) = 0, 
G = G], the above equation (|5.1|) reduces to the usual 
time evolution equation for the mass density contrast [HI, 
I771 . l80j , while in the case of the geometrical DE models 
[Q(t) = 0, G = G e ff] (eg., |zl). 

Useful expressions of the growth factor have been given 
by Peebles [80] for the ACDM cosmology. Several works 
have also derived the growth factor for w =const DE 
models f8ll - l84j ] , and for the braneworld cosmology [H[ . 
Linder [23l | and Linder and Cahn [86j derived similar ex- 
pressions for dark energy models in which the Ricci scalar 
varies with time (and which we use in the current study) , 
as well as for models with a time varying equation of 
state, while for the scalar tensor models the growth fac- 
tor was obtained by Gannouji & Polarski (87] ■ For the 
interacting DE models we use the growth factors which 
were recently derived in [lH EH. [75j. 3 



6. HALO ABUNDANCES AND THEIR 
EVOLUTION IN DARK ENERGY MODELS 

In this section we derive the cluster-size halo num- 
ber counts and their corresponding redshift distribution 
within the framework of the DE models analyzed in this 
work. We will then compare our predictions with those 
of the conventional ACDM cosmology. In principle, this 
can help us to understand better the theoretical expecta- 
tions of the current DE models, as well as to identify the 
realistic variants of the concordance ACDM cosmology. 

An important theoretical question in cosmology is how 
to determine the fraction of matter in the universe that 
has formed bounded structures, and what is its distri- 
bution in mass at any given redshift after recombination. 
The simplest successful answer to this question was given 
in 1974 by Press and Schecther [H| (hereafter PSc). In 
their approach, the primordial density fluctuations for a 
given mass M on the dark matter fluid is described by 
a random Gaussian field. The function, F(M, z), repre- 
sents the fraction of the universe that has collapsed by 
the redshift z in halos above some mass M. With this 
function one may estimate the number density of halos, 
n(M, z), with masses within the range (M, M + SM): 



n(M, z)dM 



dT(M,z) p 
dM W 



Performing the differentiation and after some algebra one 
obtains the following: 

n(M,z)dM = -±^)f PSc (a)dM 

= A^l fpSc{(j)d M, (6.2) 
M dM 1 v ' y ' 

where /psc(c) = \/2/7r(<5 c /cr) cxp(— 5 2 /2a 2 ), S c is the lin- 
early extrapolated density threshold above which struc- 
tures collapse [89j . while a 2 (M,z) is the mass variance 
of the smoothed linear density field, extrapolated to red- 
shift z at which the halos are identified. It is given in 
Fourier space by: 



a 2 (M,z) 



2tt 2 



k 2 P(k)W 2 (kR)dk , (6.3) 



where D(z) is the growth factor of perturbations, 
P{k) is the power-spectrum of the linear density field, 
W{kR) = 3(sinfci? - kRcoskR) / {kRf is the top-hat 
smoothing function which contains on average a mass 
M within a radius R = (3M / 'A-np) 1 ^ and p = 2.78 x 
10 n Sl m /i 2 Ai©Mpc~ 3 . We use the CDM power spectrum, 
P{k) = P k n T 2 (n m ,k), with T(Q m ,k) the CDM trans- 
fer function (accordingto [§3]), n — 0.96, following the 
7-year WMAP results 0], and 



P ~ 2.2 x 10" 



5 £~2 7- 



H'^k 1 



WMAP 1 



where fcwMAP = 0.02 Mpc -1 is a characteristic wave- 
length, given in [5(. Note that in this approach all the 
mass is locked inside halos, according to the normaliza- 
tion constraint: 



+ OC 



/pSc(CT)dmcr 1 = 1. 



(6.4) 



dM. 



(6.1) 



Although the above (Press-Schecther) formulation was 
shown to give a reasonable approximation to the ex- 
pectations provided by numerical simulations, it was 
later found to over-predict/under-predict the number of 
low/high mass halos at the present epoch [9ll. l92l|. 

More recently, a large number of works have provided 
better fitting functions of /(<x), some of them based on 
a phenomenological approach. In the present treatment, 
we adopt the one proposed by Reed et al. |93l |. 



3 For many DE models, it is convenient to study the growth evo- 
lution in terms of the expansion scale a or characteristic scale 
lna, rather than f . As an example, in the case of the braneworld 
cosmology, Linder & Cahn [6ll | used g = dln(D / a) / dlna (see 
their eqs. 8, 20 and 27), while for the LRDE and HRDE models, 
Linder utilized g = D/a (see eq. 31 in his paper). Since 
the pure matter universe (Einstein de-Sitter) has the solution of 
D Ed g = a, we normalize our DE models such as to get D ~ a 
at large redshifts due to the dominance of the non-rclativistic 
matter component, (e.g., 2 = 30, where the DE density is al- 
most negligible and the radiation density is less than ~ 1% of 
the matter density). 



A. Collapse Threshold and Mass Variance of DE 
models 

In order to compare the mass function predictions of 
the different DE models, it is imperative to use for each 
model the corresponding values of S c and a 2 (M,z). In 
the Press-Schecther formalism, the rms matter fluctua- 
tions is parametrized at redshift z = so that one has 
for any cosmological model: 

a(M,z)=a 8 (z) w - - } (6.5) 



8 




10" 8 
0.2 

a 

> o 

-0.2 



- - 


. 1 1 1 1 _ 

/ 


jj^^^xx- ' 

- X N 


tt ¥wrMiiil i inSii.\ 1 i " 



10 13 10 13 10 14 10 15 10 13 10 13 10 14 10 16 
M/M M/M 

FIG. 1: The halo mass function at two different redshifts. 
The different DE models are represented by different symbols 
and/or line types (see Table I for definitions). 



with 



$>(n m ,R)= k n+2 T 2 (n m ,k)W 2 (kR)dk (6.6) 
Jo 



and 



cr 8 (z) = cr 8 (0) 



D(0) ' 



■7) 



where er 8 (0)[= er 8 ] the rms mass fluctuation on i? 8 = 
Mpc scales at redshift z = 0. Therefore, for a 
given perturbed DE model, it is not enough to con- 
sider its growing mode, D(z), but one has also to ob- 
tain the corresponding erg value. Below the parame- 
ters .Dde(z), <78,de denote the growth factor and the rms 
mass perturbation normalization of the DE models, while 
D\(z),aa,A denote the corresponding quantities for the 
reference ACDM model. 



1. The relevant DE model S c values 

It is well known that in the conventional A cosmology 
S c ~ 1.675, while Weinberg & Kamionkowski [38], pro- 
vide an accurate fitting formula to estimate 5 C for any DE 
model with a constant equation of state parameter (see 
their eq.18). Let us now discuss the values of the linearly 
extrapolated density S c adopted here for the different DE 
models. 

Firstly, since for the Arc, Aps and CCDM models, the 
EoS parameter is strictly equal to —1, we are completely 
justified to use 5 C ~ 1.675. Using literature data we 
also find: (i) for the BRG model S c ~ 1.667 (ii) 



for the EDE model S c ^ 1-672 [44], (hi) for the CPL 
model, S c ~ 1.663 0, and (iv) for the VCG model, S c ~ 
1.642 [44j| . Interestingly, one can check that the above 
S c values can be well approximated using the previously 
discussed fitting formula [38| , despite the fact that it was 
derived for a constant equation of state parameter. As an 
example, in the case of the BRG geometrical model the 
fitting formula predicts S c ~ 1.667, which is in excellent 
agreement with that found by the theoretical analysis of 
Schmidt et al. 

Secondly, for the rest of the dark energy models 
(LRDE, TCM, PNGB and HRDE) there are, to our 
knowledge, no available S c values in the literature, which 
implies that one has to study in detail the spherical col- 
lapse model in order to estimate their exact S c values, 
something which is beyond the scope of the present pa- 
per. However, since in these models w ~ — 1 close to 
the present epoch, we have decided to adopt the the 
Weinberg & Kamionkowski formula because it appears 
to work quite well. In this case, for the remaining DE 
models (LRDE, TCM, PNGB and HRDE) we have de- 
rived S r ~ 1.674. 



2. Estimation of the different DE model as 

The different DE models cr 8 value can be estimated by 
scaling the present time (z = 0) ct 8 .a value to that rele- 
vant to each DE model by using eq. (|6.3p at the present 
time. We again have that: 



08, DE = 08, A' 



^De(O) P DE ,0 *(On,i? 8 



£a(o) V ^a,o *(n m , A ,i? 8 ) 



(6.8) 



with Sl m ,A denoting the value for the reference A model 
(in our case £l m ,A — 0.28, see section 3), while we have 
(for fixed Hq, as in our case), that: 



-Pde.o 
Pa.o 



1/2 



m,A 



(6.9) 



It thus follows that if the os.a value is known, one may de- 
rive the corresponding <t 8j de values to the other DE mod- 
els. In this concern, the combined SNIa+BAO-|-WMAP5 
analysis of Komatsu et al. [5] (see also 94]) pro- 
vided a value of ag,A — 0.812, while the correspond- 
ing WMAP7 analysis yielded [for w(z) = —1]: er 8 ,A — 
0.803 (using the WMAP7 alone) and ct 8 ~ 0.807 for 
the joint WMAP7+BAO+i? analysis [a]. A recent 
analysis based on cluster abundances have also fur- 
nished the following degenerate combination: <t 8i a = 
0.83 ± 0.03(On,A/0-25)-°- 41 [95|, which for our 'case 
(n m> A = 0.28) implies ct 8jA ^ 0.792. Fu et al. [H 
based on a weak-lensing procedure found ct 8i a = 0.837 ± 
0.084(ft m:A /0.25)-°- 53 which implies <t 8:A ~ 0.788 in our 
case. In addition, studies based on the peculiar ve- 
locities statistical analysis (9?J obtained ct 8i a = 0.88 ± 
0.05(ft m , A /0.25)- a53 or a s [ A ^ 0.829 in our case. It 
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should be stressed that the the average scattering of the 
four independent ct 8i a values (of those based on WMAP, 
we use only the joint WMAP7+BAO+ff result) is quite 
small (((is, a) — 0.804 ± 0.018), thereby reinforcing the 
consistency of the different measurements. Finally, by 
inserting the latter value in Eq. (|6.8[) we can estimate the 
corresponding os.de varu es, listed in Table 1, to be used 
in our mass function analysis. 

For completeness, it should be remarked that some re- 
cent analyzes are suggesting significant higher values of 
ag. For example, Watkins et al. [98| studying the bulk 
flow on scales of ~ 100 /i _1 Mpc found a erg normaliza- 
tion which is increased by a factor of ~ 2 with respect 
to the one of ACDM model. If these results are correct, 
the ACDM model would be strongly challenged. How- 
ever, such a discussion is beyond the scope of the current 
work. 



B. Halo Mass Function & Number Counts of DE 
Models 

We now pass to our results. In the upper panels of 
Fig. 1, for two different redshifts [z — 0.5 and 2), we 
display the integral halo mass function, n(> M), for all 
DE models previously discussed. The different models 
are characterized by the symbols and line- types presented 
in Table 1. 

In the corresponding lower panels one may see the frac- 
tional difference between each DE model with the concor- 
dance ACDM model, that is, Sn/n\ — (tide — n\)/nA- 
We stress that we have shown the case z = 2 only for 
comparison of the model expectations. In other words, 
it is not a statement of the viability to actually observe 
clusters at such a large redshift (not only due to tech- 
nical limitations but also because the majority of dark 
matter halos are not expected to have virialized at such 
a redshift, thereby producing their X-ray signatures). 

It is worth noticing the mass function expectation of 
the Aps model (open red circles), the BRG model (blue 
dashed line) and, at high redshifts, of the CCDM model 
(magenta stars) appear to be completely different from 
the reference ACDM cosmology. 

Given the halo mass function we can now derive an 
observable quantity which is the redshift distribution of 
clusters, Af(z), within some determined mass range, say 
Mi < M/Mq < M-2- This can be estimated by integrat- 
ing, in mass, the expected differential halo mass function, 
n(M, z), according to: 



Af(z) 



dV 

dz 



M, 



n(M, z)dM, 



(6.10) 



Mi 



where dV/dz is the comoving volume element, which in 
a flat universe takes the form: 



dV 2 dr 

-=4,r (z)-(z), 



with r(z) denoting the comoving radial distance out to 
redshift z: 



. . c f z dx 
r(z) = — / — r-. 
K > Ho J Q E{x) 



3.12) 



In Fig. 2 (upper panel), we show the theoretically ex- 
pected cluster redshift distribution, Af(z), for all the 
models studied (line and symbol types are listed in Table 
1), for cluster of galaxies size halos, ie., Mi = 10 14 M 
and M2 = 10 16 M©. It is evident that many of the 
different models show significant differences with respect 
to the concordance A model and therefore we could, in 
principle, distinguish them. In the lower panel of Fig. 2 
we show the relative differences of the various DE mod- 
els with respect to the expectations of the concordance 
ACDM model. Note that the three models (A PS , CCDM 
and BRG models), for which we found very large mass 
function differences (Fig. 1), are not shown because their 
large relative differences are beyond the limits of the 
lower panel. 

Let us now discuss the expected redshift distributions, 
based on two future cluster surveys, and also the pos- 
sibility to observationally discriminate the different DE 
models. 

These two realistic future surveys are: 

(a) the eROSITA satellite X-ray survey, with a flux limit 
of: /lim = 3.3 x 10~ 14 ergs s _1 cm -2 , at the energy band 
0.5-5 keV and covering ~ 20000 deg 2 of the sky, 

(b) the South Pole Telescope SZ survey, with a limiting 
flux density at vo = 150 GHz of f UQ ,Um — 5 mJy and a 
sky coverage of ~ 4000 deg 2 . 

To realize the predictions of the first survey we use the 
relation between halo mass and bolometric X-ra y lumi- 
nosity, as a function of redshift, provided in 100], ie: 



L{M,z) = 3.087 x 10 4 



ME(z) 
lQ 15 /i- 1 M 



1.554 



ergs 



(6.13) 

The limiting halo mass that can be observed at redshift z 
is then found by inserting in the above equation the lim- 
iting luminosity, given by: L — ^Kdj^fumCb, with dh the 
luminosity distance corresponding to the redshift z and q, 
the band correction, necessary to convert the bolometric 
luminosity of eq. (j6.13|) to the 0.5-5 keV band of eROSITA. 
We estimate this correction by assuming a Raymond- 
Smith (1977) plasma model with a metallicity of O.4Z , 
a typical cluster temperature of ~ 4 keV and a Galactic 
absorption column density of tih = 10 21 cm~ 2 . 

The predictions of the second survey can be realized 
using agai n the relation between limiting flux and halo 
mass from llOOl: 



2.592 x 10 8 mJy / M 



1.876 



.11) 



(PjjJyZ) yiu"iW0 

(6.14) 

where d^(z) = d^/ '(1 + z) 2 is the angular diameter dis- 
tance out to redshift z. 
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FIG. 3: The expected cluster redshift distribution of the dif- 
ferent DE models and for the two different future cluster sur- 
veys (upper panels), and the corresponding fractional differ- 
ence with respect to the reference ACDM model (lower pan- 
els). Symbols correspond to the different DE models as indi- 
cated in Table 1. 



differences of the DE models but only up to a redshift at 
which they are significant, that is, such that: 



A/de - A/a 



(A/DE-A/^) 1 /; 



> 3.5, 



3.15) 



FIG. 2: The expected redshift distribution of M> 1O 14 M 
clusters (upper panel) of the different DE models and the cor- 
responding fractional difference between the models and the 
reference ACDM model (lower panel). The lower panel shows 
only those DE models that have fractional relative differences, 
with respect to the ACDM model, of < 45%. Symbols corre- 
spond to the different DE models as indicated in Table 1. 



In Fig. 3 (upper panel) we present the expected red- 
shift distributions above a limiting halo mass, which is 
Mi = A/ii m it = max[10 14 M Q , Mf], with Mf correspond- 
ing to the mass related to the flux-limit at the different 
redshifts, estimated by solving eq. (|6.13j) and eq. (|6.14|) for 
M . In the lower panels we present the fractional differ- 
ence between the different DE models and ACDM, simi- 
larly to Fig. 2, but now for the realistic case of the previ- 
ously mentioned future cluster surveys. It is evident that 
the imposed flux-limits together with the scarcity of high- 
mass halos at large redshifts, induces an abrupt decline 
of Af(z) with z, especially in the case of the eROSITA X- 
ray survey (note the shallower redshifts depicted in Fig. 3 
with respect to Fig. 2). 

In the lower panels of Fig. 3 we display the relative 



with Ade the redshift distribution predicted by some 
DE model and A/a the corresponding redshift distribu- 
tion of ACDM model. However, such a criterion does 
not take into account cosmic variance and possible ob- 
servational systematic uncertainties which can hamper 
detecting small (but according to Eq. 16.151 significant) 
relative differences. We believe that relative differences 
of < 5% will be difficult to detect especially at relative 
high redshifts. 

In Table 1, one may see a more compact presentation 
of our results including the relative fractional difference 
between all DE models and the ACDM model, in 3 dis- 
tinct redshift bins and for both future surveys. 

Based on our Af(z) analysis and the results presented 
in Figure 3 and Table 1, we can now divide the studied 
DE models into those that can be distinguished obser- 
vationally and those that are practically indistinguish- 
able from the current paradigm (ACDM) . The latter DE 
models are the following four: XCDM, HRDE, TCM 
and PNGB. One has to remember however that these 
results are based on using DE model parameters that 
have been fitted by the present day cosmological data 
(see section 4). As an illustrative example, we remind 
the reader that the XCDM model, compared here, is one 
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with w — —0.99; if future cosmological data would pro- 
vide a different value for the equation of state parameter 
then the N(z) predictions of such an XCDM model could 
be quite different than those derived here. 

Regarding the models that are distinguishable with re- 
spect to the concordance model, three of them (BRG, 
Aps and CCDM) show extremely large variations mak- 
ing it trivial to distinguish them. From the rest of the 
distinguishable DE models all of them show clear signs of 
difference, at all redshifts and in both future cluster sur- 
veys, with respect to the ACDM expectations; only the 
AflG model needs to be distinguished at higher redshifts 
0^ 0.5). 

As an additional test and in order to check the sensi- 
tivity of our results only on the functional form of the DE 
equation of state parameter, we have imposed a unique 
value of the rms mass fluctuation normalization to all 
DE models according to: <j^,de = 08,A) and repeated 
the cluster-size halo N(z) analysis. We now find slightly 
different results although in the same overall direction 
with our main analysis. For example, the DE models 
that cannot be distinguished from the reference ACDM 
model are now the XCDM, TCM, PNGB, EDE and A RG , 
with the first three being common in both analyzes, the 
Aug model showing a small difference and the EDE be- 
ing the only model that shows a significantly different 
behavior between the two analyzes. 

Finally we would like to mention that an interesting 
paper appeared recently [l0l| and among other issues, it 
compares different forms of the halo mass function and 
its redshift evolution using N-body simulations of the 
ACDM and wCDM (w =const) models. They do find 
some differences at the few percent level. Although our 
analysis is self-consistent, in the sense that we compare 
the expectations of DE models with respect to those of 
the concordance cosmology using the same mass function 
model, we plan to investigate in a forthcoming paper how 
sensitive are our observational predictions to the different 
mass functions fitting formulas. 



obtaining tight constraints on the main cosmological and 
dark energy free parameters. At the level of the resulting 
Hubble function, we have found that of all dark energy 
models (apart a Brane world cosmology), are statistically 
indistinguishable (within la) from a flat ACDM model, 
as long as they are confronted with the quoted set of 
observations. 

On the other hand, despite the fact that these models 
closely reproduce the ACDM Hubble expansion, we show 
that eight out of the twelve studied DE models, using the 
observationally fitted model parameters, can be differen- 
tiated from the reference ACDM model on the basis of 
their redshift distribution of cluster-size halos. Such a 
comparison was made possible by using the expectations 
of a future X-ray (based on the eROSITA Satellite) and 
SZ cluster surveys (based on the South Pole Telescope). 

The main comparison results can be summarized in the 
following statements (for nomenclature see section 4): 

• Four DE models, namely, XCDM, HRDE, TCM 
and PNGB cannot be distinguished from the 
ACDM model at any significant level. 

• Seven models, ie., the BRG, CPL, LRDE, EDE, 
VCG, Aps, CCDM, can be easily distinguished due 
to the fact that they show strong and significant 
variations with respect to the concordance A model 
even at z = 0, implying that even with the present 
day surveys one could effectively distinguish them. 

• The Ape model, although presents relatively small 
variations with respect to the concordance A model, 
it can be clearly distinguished at relatively high 
redshifts (z^ 0.5). 

In a future work we will present a comparison be- 
tween our model predictions and the observationally de- 
termined cluster mass function at different redshifts as 
well as the available (X-ray or optical) cluster redshift 
distribution. 



7. CONCLUDING REMARKS 

In this paper, we have investigated the cluster abun- 
dances beyond the conventional ACDM cosmology by us- 
ing several parameterizations for the dark energy. In 
order to do that, we first performed a joint likelihood 
analysis using the most recent high quality cosmological 
data (SNIa, CMB shift parameter and BAOs), thereby 
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TABLE I: Numerical results. The 1 st column indicates the DE model. 2 nd column are symbols or line types of the models 
appearing in Figs. 2, 3 and 4. 3 rd and 4 th columns show the <T8,de and 8 C values. The remaining columns present the fractional 
relative difference between the models and the ACDM cosmology for two future cluster surveys discussed in the text. Bold 
letters denote those models that can be clearly distinguished from fiat ACDM model at some redshift bin (uncertainties appear 
only when their value is > 10~ 2 ). 
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The mass function of cluster-size halos and their redshift distribution are computed for 12 distinct 
accelerating cosmological scenarios and confronted to the predictions of the conventional flat ACDM 
model. The comparison with ACDM is performed by a two-step process. Firstly, we determine the 
free parameters of all models through a joint analysis involving the latest cosmological data, using 
SNe type la, the CMB shift parameter and BAO. Apart from a brane world inspired cosmology, 
it is found that the derived Hubble relation of the remaining models reproduce the ACDM results 
approximately with the same degree of statistical confidence. Secondly, in order to attempt distin- 
guish the different dark energy models from the expectations of ACDM, we analyze the predicted 
cluster-size halo redshift distribution on the basis of two future cluster surveys: (i) an X-ray sur- 
vey based on the eROSITA satellite, and (ii) a Sunayev-Zeldovich survey based on the South Pole 
Telescope. As a result, we find that the predictions of 8 out of 12 dark energy models can be 
clearly distinguished from the ACDM cosmology, while the predictions of 4 models are statistically 
equivalent to those of the ACDM model, as far as the expected cluster mass function and redshift 
distribution are concerned. The present analysis suggest that such a technique appears to be very 
competitive to independent tests probing the late time evolution of the Universe and the associated 
dark energy effects. 
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1. INTRODUCTION 

Recent studies in observational cosmology lead to the 
conclusion that the available high quality cosmological 
data, from Supernovae type la (SNe la), matter power 
spectrum analysis, the angular power spectrum of the 
cosmic microwave background (CMB), baryon acoustic 
oscillations (BAO) and other complementary probes, like 
the existence of old galaxies at high rcdshifts, are well fit- 
ted by an emerging "standard cosmological model" (see 
[1-8] and references therein). By assuming flatness, as 
theoretically required by inflation and the CMB observa- 
tions, the standard model is described by the Friedmann 
equation: 



H 2 (a) = (- J = — [p m (a) + PDE (a)} , (1.1) 

where a(t) is the scale factor of the universe, p m {a) is the 
density corresponding to the sum of baryonic plus cold 
dark matter (the latter needed to explain clustering) , and 
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an extra component Pbe((i), with negative pressure called 
dark energy (DE hereafter), that accelerates the cosmic 
expansion (for reviews see [9]). 

Nowadays, the physics of dark energy is considered one 
of the most fundamental and challenging problems on 
the interface uniting Astronomy, Cosmology and Parti- 
cle Physics. In the last decade there have been theoret- 
ical debates among the cosmologists regarding the na- 
ture of this exotic component. Many candidates have 
been proposed in the literature, such as a cosmolog- 
ical constant A (vacuum), time- varying A(t) cosmolo- 
gies, quintessence, k— essence, quartessence, vector fields, 
phantom, tachyons, modifications of gravity, Chaply- 
gin gas and the list goes on (see [10-26] and references 
therein) . Gcncrically, some high energy field theories also 
suggest that the equation of state of such a dark energy 
is a function of cosmic time (see, for instance, [27]). Nat- 
urally, in order to establish the evolution of its equation 
of state (EoS), a realistic form of H(a) is required and 
must be constrained through a combination of indepen- 
dent DE probes. 

On the other hand, the abundance of collapsed struc- 
tures as a function of mass and redshift is a key statisti- 
cal test for studies of matter distribution in the universe, 
and, more importantly, it can readily be accessed from 
observations [28]. Indeed, the mass function of galaxy 
clusters has been measured based on X-ray surveys [29- 



31], via weak and strong lensing studies [32-34], using op- 
tical surveys, like the SDSS [35, 36], as well as, through 
Sunayev-Zeldovich (SZ) effect [37]. In the last decade 
many authors have been involved in this kind of stud- 
ies and have found that the abundance of the collapsed 
structures is affected by the presence of a dark energy 
component [38-46]. 

In this work, we discuss how to differentiate a large 
family of flat DE cosmologies (12 models) from the con- 
ventional ACDM model. Initially, a joint statistical anal- 
ysis, involving the latest observational data from SNe 
type la, CMB shift parameter and BAO, is implemented. 
Since however the resulting Hubble functions of these DE 
models (apart a brane world scenario) are statistically in- 
distinguishable, we attempt to discriminate the different 
DE models by computing the halo mass function and the 
corresponding redshift distribution of the cluster-size ha- 
los. Finally, by using future X-ray and SZ surveys we 
show that the evolution of the cluster abundances (espe- 
cially at large redshifts 0.5), is a potential discrimina- 
tor (in the Popper sense) for a large fraction of the stud- 
ied DE models. As an extra bonus, we find that many 
of the DE models could be differentiated, from ACDM, 
even with the present available observational data from 
galaxy cluster mass function. 

The article is planned as follows. Section 2 is a pro- 
vision of dark energy issue and the overall approach 
adopted in the paper. In section 3, a joint statistical 
analysis based on SNe la, CMB and BAO is carried out 
for the conventional ACDM model. This statistical ap- 
proach is used to constrain the DE model parameters and 
it is presented in section 4. The linear growth factor of all 
models is discussed in section 5, while in section 6, we dis- 
cuss and compare the corresponding theoretical predic- 
tions regarding the evolution of the cluster abundances. 
Finally, the main conclusions are summarized in section 
7. Throughout the paper we adopt H = 71km/sec/Mpc. 



EoS parameter in terms of E(a) = H(a)/H [47] using 
the Friedmann equations as 



w(a) 



-1 



2 tflnE 
3 U da 



1 - Q m a- 3 E~ 2 (a)' 



(2.2) 



Since the exact nature of the DE has yet to be found, 
the above DE EoS parameter encodes our ignorance re- 
garding the physical mechanism powering the late time 
cosmic acceleration. 

The methodology described above can also be applied 
to the framework of modified gravity (see [23, 48]). In 
this case, instead of using the exact Hubble flow through 
a modification of the Friedmann equation one may con- 
sider an equivalent Hubble flow somewhat mimicking Eq. 
(1.1). The key point here is that the accelerating ex- 
pansion can be attributed to a kind of "geometrical" 
DE contribution. Now, since the matter density (bary- 
onic+dark) cannot accelerate the cosmic expansion, we 
perform the following parametrization [23, 48]: 



E 2 (a) 



H 2 (a) 
til 



SH 2 



(2.3) 



Naturally, any modification to the Friedmann equation 
of general relativity may be included in the last term of 
the above expression. After some algebra one may also 
derive, using Eqs. (2.2) and (2.3), an effective ("geomet- 
rical") dark energy EoS parameter, given by: 



w(a) 



1 dlnSH 2 
3 rflna 



(2.4) 



The above formulation will be adopted in our statistical 
analysis of all DE models discussed in section 4. 



3. LIKELIHOOD ANALYSIS: THE ACDM CASE 



2. DARK ENERGY EQUATION OF STATE 

In what follows, it will be assumed that the universe 
is a self-gravitating fluid described by general relativ- 
ity, and endowed with a spatially flat homogeneous and 
isotropic geometry. In addition, we also consider that 
it is filled by non-relativistic matter plus a DE compo- 
nent (or some effective mechanism that simulates it), and 
whose equation of state, pde = w(a)pDE, is driving the 
present accelerating stage. Following standard lines, the 
Hubble parameter reads: 



(2.1) 



where E(a) is the normalized Hubble flow, TO is the di- 
mensionless matter density at the present epoch, Dde = 
1 — Cl m denotes the DE density parameter, and w(a) its 
EoS parameter. On the other hand, we can express the 



In this section wc briefly present the basic observa- 
tional samples and statistical analysis tools that will be 
used to constrain the cosmological parameters of the DE 
models. Here we discuss the ACDM model since it is now 
widely believed that if a better cosmology is called for, 
it will describe a cosmos that looks much like a ACDM 
model. 

(i) Supernovas. In our statistical analysis we consider 
the Constitution set, containing 397 SNe type la, as com- 
piled by Hicken et al. [6]. In order to avoid possible 
problems related to the local bulk flow, we use a subset 
of this sample containing 366 SNe la all with redshifts 
z > 0.02. The likelihood estimator is determined by a 
XsNia statistics 



366 r 



XSNIa(P) = 



i=l 



M th K,p) - M obs (ai) 



(3.1) 



where dj = (1 + z{) 1 is the scale factor of the Universe 
in the observed redshift Zi, fi is the distance modulus 



fj, = m — M = 51ogrfi + 25 and (a, p) is the luminosity 
distance, d L (a,p) = ca- 1 p^y. 

(ii) Baryon Acoustic Oscillations (BAO). In addition 
to the SNe la data, we also consider the BAO scale 
produced in the last scattering surface by the competi- 
tion between the pressure of the coupled baryon-photon 
fluid and gravity. The resulting acoustic waves leaves 
(in the course of the evolution) an overdensity signature 
at certain length scales of the matter distribution. Ev- 
idence of this excess was recently found in the cluster- 
ing properties of the SDSS galaxies [49], and it provides 
a suitable "standard ruler" for constraining DE mod- 
els. In particular, we consider the following estimator 

^(P) = [z*E(a^)]i/3 fa, 



2/3 



measured from the 



SDSS data to be A = 0.469 ± 0.017, where z s = 0.35 [or 
a s = (1 + Zg)^ 1 ~ 0.75]. Therefore, the corresponding 
Xbao function can be written as 



Xbao(p) 



[A(p) - 0.469] 2 
0.017 2 



(3.2) 



where p is a vector containing the cosmological fitting 
parameters. 

(iii) CMB Shift Parameter. Another interesting geo- 
metrical probe for dark energy is provided by the an- 
gular scale of the sound horizon at the last scattering 
surface. It is encoded in the location l\ T of the first 
peak of the angular (CMB) power spectrum [50, 51], and 
may be defined by the quantity 1Z = y/£l m a -i d ^ a ^ ■ 
The shift parameter measured from the WMAP 7-years 
data [5] is K = 1.726 ± 0.019 at z Ls = 1091.36 [or 
ai s = (1 + zis)^ 1 ~ 9.154 x 10~ 4 ]. In this case, the 
Xcmb function reads 



Xcmb(P) 



[R(p) - 1.726] 2 



(3.3) 



0.019 2 

It should be stressed that for CMB shift parameter, the 
contribution of the radiative component, (f2#a -4 , where 
Q R ~ 4.174 x 10 _5 /i -2 ) needs also to be considered [5]. 
Note also that the measured CMB shift parameter is 
somewhat model dependent but mostly to models which 
arc not included in our analysis. For example, such is the 
case when massive neutrinos are included or when there 
is a strongly time varying equation of state parameter. 
The robustness of the shift parameter has been tested 
and discussed in [52]. 

In order to put tighter constraints on the correspond- 
ing parameter space of each DE model, the probes de- 
scribed above must be combined through a joint like- 
lihood analysis 1 , given by the product of the individual 



likelihoods according to: £ to t(p) = £sNia x ^bao x £ cm b, 
which translates in the joint x 2 function in an addition: 
xL(p) = XsNia + Xbao + Xcmb- The results of this anal- 
ysis provide: TO = 0.28 ± 0.01 with Xtot(^m)/dof ~ 
439.6/367. 

In this concern, it should be remarked that a value of 
XsNia — 439 is what one gets by using directly the Con- 
stitution SNIa set of Hicken et al. [6] since the number of 
the SNIa used (366) dominates the overall Xtot budget of 
the joint likelihood. In addition, the corresponding good- 
ness of fit (XsNia/dof — 1-2) is significantly larger then 
the one found by Davis et al. [2] (xsNia/^of — 1) from a 
previous sample containing 192 SNIa (see the discussion 
in [53, 54]). Such a discrepancy appears to be the out- 
come of the different approaches chosen in order to join 
the different contributing SNIa sets. According to Hicken 
(private communication, 2009) in the case of the Davis 
et al. [2] data the nearby SNIa were imposed to provide 
X 2 /dof ~ 1 by hand, while no such fine-tuning was im- 
posed on the Constitution set. Also the latter set includes 
distant SNIa which have typically larger distance modu- 
lus uncertainties with respect to those used in Davis et 
al. sample [2]. In particular, this means that the higher 
XsNia/dof value of the Constitution set should probably 
be attributed to typically lower SNIa distance modulus 
uncertainties. 2 

By using the most recent BAO results of Percival et 
al. [55], we have checked and verified that the above 
constraints of the flat ACDM model are not apprecia- 
bly modified. Therefore, the reader should keep in mind 
the value Xtot — 439.6 for any further comparison with 
results predicted by the alternative DE models. 



4. CONSTRAINTS ON DARK ENERGY 
MODELS 

Let us now present the twelve flat DE models whose 
free parameters will be constrained by using the same 
methodology and cosmological data as that applied to 
the ACDM model (see previous section). 



A. Constant Equation of State (XCDM model) 

In this kind of cosmology (hereafter XCDM-modcls) 
the equation of state parameter is constant [56]. In a 
point of fact, these DE models do not have a clear phys- 
ical motivation. In particular, for quintessence models 
driven by a real scalar field, a constant EoS parameter 
requires an extreme fine tuning of its potential and ki- 
netic energy. In spite of that, this subclass of DE models 



1 Likelihoods are normalized to their maximum values. In the 
present analysis we always report la uncertainties on the fitted 
parameters. Note also that the total number of data points used 
here is Ntot = 368, while the associated degrees of freedom are: 
dof= Ntot — "fin with rtfl t the model-dependent number of fitted 
parameters. 



2 Such a possibility has been crudely tested by Plionis et al. [54]. A 
lower value Xs^j a /dof ~ 1.07 (similar to that of Davies [2]) has 
been obtained by increasing the distance modulus uncertainty of 
the Constitution nearby SNIa (z < 0.4) by 20%. 



have been widely used in the literature due to their sim- 
plicity. Notice that DE models with a canonical kinetic 
term have — 1 < w < —1/3. Models with (w < —1), 
sometimes called phantom dark energy [57], are endowed 
with a very exotic nature, like a scalar field with negative 
kinetic energy. Now, by using Eq. (2.1) the normalized 
Hubble parameter becomes 

E 2 (a) = n m a-* + (1 - n m )a-^ 1+w \ (4.1) 

Now, in order to constrain XCDM-models with the 
observational data we sample Q m G [0.1,1] and w G 
[—2, —0.4] in steps of 0.01. As a result, we find that the 
best fit values are Q m = 0.28±0.01 and w = -0.99±0.05 
with Xtot( n m,w)/dof ~ 439.5/366. 

The above results are in excellent agreement with those 
found by different authors [2, 4-6, 8, 31]. It is also worth 
noticing that the concordance ACDM cosmology is de- 
scribed by a XCDM model with w strictly equal to -1. 
The corresponding limits form the basis of our present 
comparison and were separately derived (see last section) . 



B. Parametric Dark Energy (CPL model) 

This kind of cosmology was first discussed with basis 
on the Chcvalicr-Polarski-Lindcr [58, 59] paramctrization 
(hereafter CPL). The dark energy EoS parameter is de- 
fined as a first order Taylor expansion around the present 
epoch: 



w(a) = wq + u>i(l — a), 



(4.2) 



where Wo and w± are constants. The normalized Hubble 
parameter now reads: 



E 2 (a) = n m a 3 + (l-fi TO )a 



-3(l+wo+«'i) (a 3u)i(a— 1) 



(4.3) 



In order to constrain the free parameters, we sample 
them as follows: wo € [—2, —0.4] and w\ G [—2.6, 2.6] 
in steps of 0.01. By fixing a prior in the density parame- 
ter, Vt m = 0.28, we find that the overall likelihood func- 
tion peaks at w a = -0.96 ± 0.13 and w x = -0.40 ± 0.70 
in very good agreement with the 7 years WMAP data 
[5]. The corresponding Xt o t( w 0,wi)/dof is 439.5/366. A 
value that should be compared with the ACDM predic- 
tion. 



functional form of the normalized Hubble flow as given 
by Eq. (2.3) is affected only by one free parameter, fl m . 
The quantity SH 2 is given by 



6H 2 



(4.4) 



where flbw = (1 - ^m) 2 /4. From Eq.2.4, it is readily 
checked that the geometrical DE equation of state pa- 
rameter reduces to 



w(a) — 



1 



1 + fi m (a) 



(4.5) 



where fi m (a) = tt m a~ 3 /E 2 (a). The joint likelihood anal- 
ysis provides a best fit value of TO = 0.23 ±0.01, but the 
fit is much worse, Xtot(^m)/dof ~ 500/367, in compari- 
son with the one provided by a ACDM cosmology. 



D. Low Ricci Dark Energy (LRDE model) 

In this model we use a simple parametrization for the 
Ricci scalar which is based on a Taylor expansion around 
the present time: 1Z(a) = r +ri(l — a). It is interesting 
to mention that at the early epochs the cosmic evolution 
tends asymptotically to be matter dominated (for more 
details see [23]). In this framework, we have 



E 2 (a) 



a 4(r +ri-l) e 4ri(l-o) 



a> a t 
a < a t 



(4.6) 



where at = 1 — (1 — 4ro)/4ri. At present, the matter 
density parameter is directly related with the above con- 
stants via 



4r + 4r : 



1 



4n 



4r +4ri- 



3 l-4r 



The corresponding EoS parameter is given by 



w(a) 



1 - 411(a) 



i-q.-/:^ 1 -^) 



(4.7) 



(4.8) 



Notice, that we sample the unknown parameters as fol- 
lows: r G [0.5,1.5] and n G [—2.4,-0.1] and here are 
the results: r = 0.80 ± 0.02 and n = -0.69 ± 0.05 
(Q m ~ 0.27) with xlt(ro,ri)/dof ~ 439.9/366. 



E. High Ricci Dark Energy (HRDE model) 



C. Braneworld Gravity (BRG model) 

In the context of a braneworld cosmology (hereafter 
BRG) the accelerated expansion of the universe can be 
explained by a modification of the gravitational interac- 
tion in which gravity itself becomes weak at very large 
distances (close to the Hubble scale) due to the fact that 
our four dimensional brane survives into an extra dimen- 
sional manifold (see [60] and references therein). The in- 
teresting aspect of this scenario is that the corresponding 



A different Ricci DE model was proposed by Linder 
[61]. In this framework, the Ricci scalar evolves at high 
redshifts obeying the following expression 



1 — 3wt 



SH 2 



H 2 



(4.9) 



where SH 2 — E 2 (a) — Vt m a 3 . In this model the normal- 
ized Hubble flow becomes: 

_ n „-3 d i fl„-3u> ^- lnS "Wln(l+/3) 



E 2 (a) = n m a- 3 (1 + [3a- 3w °y 



(4.10) 



where f3 = ri" 1 — 1. As in the previous case, the effec- 
tive equation of state parameter is again related to E(a) 
according to Eq.(2.2). Now, by using the same sampling 
as in the XCDM-modcls, we find that the joint likelihood 
peaks at O m = 0.28 ± 0.01 and w a = -0.99 ± 0.05 with 
w ) I 'do / ~ 439.5/366. 



F. Tension of Cosmic Magnetic Fields 
(TCM model) 

A couple of years ago, Contoupolos and Basilakos [62] 
proposed a novel idea which is based on the following con- 
sideration (hereafter TCM): if the cosmic magnetic field 
is generated in sources (such as galaxy clusters) whose 
overall dimensions remain unchanged during the expan- 
sion of the Universe, the stretching of this field by the 
cosmic expansion generates a tension (negative pressure) 
that could possibly explain a small fraction of the DE 
(<~ 2 — 5%). In this flat cosmological scenario the nor- 
malized Hubble flow becomes: 

E 2 (a) = n m a- 3 + SH 2 , SH 2 = n A + n B a- 3+n (4.11) 

where f2s is the density parameter for the cosmic mag- 
netic fields and f^A = 1 — — ^b- From Eq. (2.4) it 
is readily seen that the effective EoS parameter which is 
related to magnetic tension reads 



w(a) = — 



3J7a + nQBd 

3(f^A + ^Bfl' 



-3+n 
-3+n\ 



(4.12) 



Again, in order to constrain the parameters, we sample 
Qb € [0.,0.3] and n e [0,10]. By considering a prior of 
n m = 0.28, the best fit values are: fl B = 0.01 ± 0.005 
and n = 1.80 ± 0.80 with Xtot(^B,n)/dof ~ 439.8/366. 



G. Pseudo-Nambu Goldstone Boson 
(PNGB model) 

In this cosmological model the DE equation of state 
parameter is expressed with the aid of the potential 
V{4>) oc [1 + «m#/F)] [63]: 



w(a) 



-1 + (I + Wq)^ 



(4.13) 



where w G [-2,-0.4] and F <E [0,8]. We see that for 
a«l, the EoS parameter goes to w(a) = —1. Based on 
this parametrization the quantity SH 2 (see eq.2.4) takes 
the following form 



SH 2 = (1 - Q m )exp 



3(l + «J2) (1 _ a , ) 



(4.14) 



In this case, our joint statistical analysis yields that the 
likelihood function peaks at wo = — 0.97± 0.09 and F = 
5.9 ± 3.2 with xtot( w o, F)/dof ~ 439.7/366. 



H. Early Dark Energy (EDE model) 

Another cosmological scenario that we include in our 
paper is the early dark energy model (hereafter EDE). 
The basic assumption here is that at early epochs the 
amount of DE is not negligible [64]. In this framework, 
the overall DE component is given by 



^de(o) 



1 



fie(l 



-3id 



i - n m + n m a 3w ° 

(4.15) 

where fl e is the early DE density and wo is the equation 
of state parameter at the present epoch. We observe that 
the EDE model was designed to satisfy the two following 
conditions: (i) mimic the effects of a late dark energy, 
and (ii) provide a decelerating expansion of the Universe 
at high redshifts. The normalized Hubble parameter is 
now written as 



E 2 (a) 



a 



1 - ^DE(a) 



(4.16) 



and from Eq. (2.2), one may obtain the EoS parameter 
as a function of the scale factor. 

Now, from the joint likelihood analysis we find that 



0.02 ± 0.02 and w 

2 



-1.04 



+0.07 



n (for the prior 

of n m = 0.28) with xtot( n e,wo)/dof ~ 439.3/366. Note, 
that a different class of EDE models was recently studied 
in [45]. 



I. Variable Chaplygin Gas (VCG model) 

Accelerating cosmologies can also be driven by the so 
called variable Chaplygin gas (hereafter VCG) which cor- 
responds to a Born- Infold tachyon action [21, 65]. Re- 
cently, an interesting family of Chaplygin gas models 
was found to be consistent with the current observational 
data [66] . In the framework of a spatially flat geometry, it 
can be shown that the normalized Hubble function takes 
the following formula: 



E 2 (a) = n b a- 3 +(l-n b )^B s a- e + (1 - B s )a~ n (4.17) 

where ~ 0.0226/i -2 is the density parameter for the 
baryonic matter [5] and B s e [0.01, 0.51] and n G [-4, 4]. 
The effective matter density parameter is: tt m = Oj, + 
(1 — VLb)^/B~ s . We find that the best fit parameters are 
B s = 0.05 ± 0.01 and n = 1.65±g ^§ (fi m ~ 0.26) with 
xl t (B s ,n)/dof ~ 441.3/366. 



J. Time Varying Vacuum (Arc and Aps models) 

Let us now consider the possibility of a decaying A- 
cosmology, that is, A = A(a). In this kind of model, the 
value of A is small because the Universe is too old, and, 
therefore, it alleviates the cosmological constant prob- 
lem. As it appears, the decaying vacuum equation of 



state does not depend on whether A is strictly constant 
or variable. Therefore, the EoS takes the usual form, 
P A (t) = -p\(t) = -A(t)/8nG. The global dynamics of 
such models have been investigated extensively in the lit- 
erature, in fact, much before the discovery of the present 
accelerating stage [13] (for a recent and extensive dis- 
cussion see Basilakos et al. [14]). By introducing in the 
Friedman equations the idea of a time-dependent vacuum 
one obtains 



A 
2 



(4.18) 



The traditional A = const cosmology can be described 
directly by the integration of Eq. (4.18), but the same 
equation is also valid for A = A (a). This means that 
a supplementary equation for the time evolution of A is 
needed in order to unveil the dynamics of this pattern. In 
this work we consider two different versions of the A (a) 
models, namely renormalization group of the quantum 
field vacuum [25, 67] (hereafter Ape) and a power series 
vacuum [68] (hereafter Aps) 



A(H) = 



A + 3 7 (H 2 - H 2 ) 
Hoh - 3Q m )H + (3 - 7 )# 2 



for A RG 
for Aps 

(4.19) 

Naturally, the vacuum energy density is normalized to its 
present value: A = A(_ff ) = 3Q\Hq (Q\ = 1 — fl m ). 
Inserting the above expressions into Eq. (4.18) we finally 
obtain the normalized Hubble flow 



E(a) = 



i-fi„ 



l- 7 
1 



30^ 

7 



-a 

7 

7 



3(1-7) 



1/2 



for Ap G 
for Aps 

(4.20) 

It has been shown [67], [70] that for A(t) there is a cou- 
pling between the time-dependent vacuum and matter 
component. Indeed, by combining the conservation of 
the total energy with the variation of the vacuum en- 
ergy, one can shown that the A(t) models provide either 
a particle production process or that the mass of the dark 
matter particles increases [69]. 

Now, by applying our statistical procedure for both 
cosmologies, the best fit parameters are: 

Ap G model: fl m = 0.28 ± 0.01 and 7 = 0.002 ± 0.001 
with x 2 tot (n mr/ )/dof ~ 439.5/366. 

Ap S model: fl m = 0.32 ± 0.02 and 7 = 3.45^^ with 
X L(^,7)M>/- 440.7/366. 



K. Creation Cold Dark Matter (CCDM Model) 



dark matter particles), and the energy conservation laws 
for each component have been investigated thoroughly by 
[72-74]. In this framework, the basic equation governing 
the global dynamics of a Universe endowed with a flat 
geometry is [72-75] 



H+-H* 



^Gp dm T_ 
3 H' 



(4.21) 



where pd m is the dark matter energy densities and T is 
the so called creation rate of the cold dark matter and 
it has units of (time) -1 . The creation pressure is nega- 
tive and defined in terms of the creation rate and other 
cosmological quantities by the expression [72-75] 



Pc = — 



PdmF_ 

3H ' 



(4.22) 



In CCDM models, the functional form of T is phcnomcno- 
logically parametrized [74] by the following relation 



3 n A (p^l) h , 

\PdmJ 



(4.23) 



where (l\ (called a in the [74]) is a constant, p co = 
3Hq/8ttG is the present day value of the critical den- 
sity. Performing now the integration of Eq. (4.21) we 
obtain 



E 2 (a) = £l m a~ 3 + Qj, 



(4.24) 



where &a = 1 — i\n- Interestingly, the CCDM model re- 
sembles the global dynamics of the concordance ACDM 
cosmology by including only one free parameter. In 
particular, this means that TO = 0.28 ± 0.01 with 
X 2 tot (n m )/dof ~ 439.6/367. 



5. EVOLUTION OF MATTER 
PERTURBATIONS 

In this section we generalize the basic equations which 
govern the behavior of the matter perturbations within 
the framework of the previously described DE models. 
We focus our attention to the most generic models, ie., 
the A(t) and CCDM models. In these cases, a neo- 
Newtonian description can be introduced by considering 
an extended continuity equation together with the Euler 
and Poisson equations [76, 79]. In virtue of the parti- 
cle creation model, the evolution equation of the growth 
factor becomes [14, 75]: 



In the context of the standard general relativity the- 
ory, is also possible to reduce the dark sector by consid- 
ering the presence of the gravitationally induced parti- 
cle creation mechanism. Like A(t)-models this kind of 
scenario has also been discussed long ago [71]. More re- 
cently, the cosmological equations for the mixture of ra- 
diation, baryons and cold dark matter (with creation of 



D + (2H + Q)D - (ATtGp m - 2HQ - q) D = 0, (5.1) 

where p m is the matter density, G — G and 

n(A _ / -A/pm for A jig or A P5 (t . 

WJ 1 3n APco H/p m for CCDM 



In the case of non interacting DE models, [Q(t) = 0, 
G = G], the above equation (5.1) reduces to the usual 
time evolution equation for the mass density contrast [48, 
77, 80], while in the case of the geometrical DE models 
[Q(t) = 0, G = G cff ] (eg., [78]). 

Useful expressions of the growth factor have been given 
by Peebles [80] for the ACDM cosmology. Several works 
have also derived the growth factor for w =const DE 
models [81-84], and for the braneworld cosmology [85]. 
Lindcr [23] and Lindcr and Cahn [86] derived similar ex- 
pressions for dark energy models in which the Ricci scalar 
varies with time (and which we use in the current study), 
as well as for models with a time varying equation of 
state, while for the scalar tensor models the growth fac- 
tor was obtained by Gannouji & Polarski [87]. For the 
interacting DE models we use the growth factors which 
were recently derived in [14, 41, 75]. 3 



6. HALO ABUNDANCES AND THEIR 
EVOLUTION IN DARK ENERGY MODELS 

In this section we derive the cluster-size halo num- 
ber counts and their corresponding rcdshift distribution 
within the framework of the DE models analyzed in this 
work. We will then compare our predictions with those 
of the conventional ACDM cosmology. In principle, this 
can help us to understand better the theoretical expecta- 
tions of the current DE models, as well as to identify the 
realistic variants of the concordance ACDM cosmology. 

An important theoretical question in cosmology is how 
to determine the fraction of matter in the universe that 
has formed bounded structures, and what is its distri- 
bution in mass at any given redshift after recombination. 
The simplest successful answer to this question was given 
in 1974 by Press and Schccther [88] (hereafter PSc). In 
their approach, the primordial density fluctuations for a 
given mass M on the dark matter fluid is described by 
a random Gaussian field. The function, F(M,z), repre- 
sents the fraction of the universe that has collapsed by 
the redshift z in halos above some mass M. With this 
function one may estimate the number density of halos, 
n(M, z), with masses within the range (M, M + SM): 



n(M, z)dM = 



df{M,z) p 
dM M 



dM. 



(6.1) 



Performing the differentiation and after some algebra one 
obtains the following: 



n(M, z)dM 



p / 1 da\ 
~M \adMJ 
p diner _1 



M dM 



f PSc (a)dM 
f PSc (a)dM, (6.2) 



where /psc(c) = \J2/-k(8 c / a) exp(— 5 2 /2a 2 ), S c is the lin- 
early extrapolated density threshold above which struc- 
tures collapse [89] , while a 2 (M, z) is the mass variance 
of the smoothed linear density field, extrapolated to red- 
shift z at which the halos are identified. It is given in 
Fourier space by: 



CT 2 (M,z) 



D 2 (z) 
2tt 2 



k 2 P{k)W 2 {kR)dk 



(6.3) 



where D(z) is the growth factor of perturbations, 
P(k) is the power-spectrum of the linear density field, 
W{kR) = 3(sinfci? - kRcoskR) / (kRf is the top-hat 
smoothing function which contains on average a mass 
M within a radius R = (3M/47rp) 1/3 and p = 2.78 x 
lO n £! m /i 2 M0Mpc~ 3 . We use the CDM power spectrum, 
P{k) = P k n T 2 (n m ,k), with T(Q m ,k) the CDM trans- 
fer function (according to [90]), n ~ 0.96, following the 
7-year WMAP results [5], and 



P ~ 2.2 x 10" 



5 



TT-Ai\-n 
"0 ''WMAP J 



where /cwmap = 0.02 Mpc -1 is a characteristic wave- 
length, given in [5]. Note that in this approach all the 
mass is locked inside halos, according to the normaliza- 
tion constraint: 



/+oo 
-oo 



/pSc(o-) dhicr 1 = 1. 



(6.4) 



Although the above (Press-Schecther) formulation was 
shown to give a reasonable approximation to the ex- 
pectations provided by numerical simulations, it was 
later found to over-predict/under-predict the number of 
low/high mass halos at the present epoch [91, 92]. 

More recently, a large number of works have provided 
better fitting functions of /(cr), some of them based on 
a phenomenological approach. In the present treatment, 
we adopt the one proposed by Reed et al. [93]. 



3 For many DE models, it is convenient to study the growth evo- 
lution in terms of the expansion scale a or characteristic scale 
lna, rather than t. As an example, in the case of the braneworld 
cosmology, Linder & Cahn [61] used g = dln(D / a) / dlna (see 
their eqs. 8, 20 and 27), while for the LRDE and HRDE models, 
Linder [86] utilized g = D/a (see eq. 31 in his paper). Since 
the pure matter universe (Einstein de-Sitter) has the solution of 
£*EdS = a > we normalize our DE models such as to get D ~ a 
at large redshifts due to the dominance of the non-rclativistic 
matter component, (e.g., z = 30, where the DE density is al- 
most negligible and the radiation density is less than ~ 1% of 
the matter density). 



A. Collapse Threshold and Mass Variance of DE 
models 

In order to compare the mass function predictions of 
the different DE models, it is imperative to use for each 
model the corresponding values of 5 C and cr 2 (M, z). In 
the Press-Schecther formalism, the rms matter fluctua- 
tions is parametrized at redshift z = so that one has 
for any cosmological model: 

a(M,z)=a 8 (z) w - - } (6.5) 
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FIG. 1: The halo mass function at two different redshifts. 
The different DE models are represented by different symbols 
and/or line types (see Table I for definitions). 



with 



*(n m ,-i) 



k n+2 T 2 (Q m ,k)W 2 (kR)dk (6.6) 



and 



a 8 (z) = er 8 (0) 



D(z) 
D(0) ' 



(6.7) 



where <7g(0) [= erg] the rms mass fluctuation on i?s = 
8/i _1 Mpc scales at redshift z = 0. Therefore, for a 
given perturbed DE model, it is not enough to con- 
sider its growing mode, D(z), but one has also to ob- 
tain the corresponding ag value. Below the parame- 
ters -Dde(-z), CT8,de denote the growth factor and the rms 
mass perturbation normalization of the DE models, while 
Da(z), o^a denote the corresponding quantities for the 
reference ACDM model. 



1. The relevant DE model 8 C values 

It is well known that in the conventional A cosmology 
S c ~ 1.675, while Weinberg & Kamionkowski [38], pro- 
vide an accurate fitting formula to estimate S c for any DE 
model with a constant equation of state parameter (see 
their eq.18). Let us now discuss the values of the linearly 
extrapolated density S c adopted here for the different DE 
models. 

Firstly, since for the Arc, Aps and CCDM models, the 
EoS parameter is strictly equal to —1, we are completely 
justified to use S c ~ 1.675. Using literature data we 
also find: (i) for the BRG model S c ~ 1.667 [99], (ii) 



for the EDE model, S c ~ 1.672 [44], (iii) for the CPL 
model, S c ~ 1.663 [44], and (iv) for the VCG model, S c ~ 
1.642 [44]. Interestingly, one can check that the above 
6 C values can be well approximated using the previously 
discussed fitting formula [38] , despite the fact that it was 
derived for a constant equation of state parameter. As an 
example, in the case of the BRG geometrical model the 
fitting formula predicts S c ~ 1.667, which is in excellent 
agreement with that found by the theoretical analysis of 
Schmidt et al. [99]. 

Secondly, for the rest of the dark energy models 
(LRDE, TCM, PNGB and HRDE) there are, to our 
knowledge, no available S c values in the literature, which 
implies that one has to study in detail the spherical col- 
lapse model in order to estimate their exact 5 C values, 
something which is beyond the scope of the present pa- 
per. However, since in these models w ~ — 1 close to 
the present epoch, we have decided to adopt the the 
Weinberg & Kamionkowski formula because it appears 
to work quite well. In this case, for the remaining DE 
models (LRDE, TCM, PNGB and HRDE) we have de- 
rived S c ~ 1.674. 



2. Estimation of the different DE model as 

The different DE models a% value can be estimated by 
scaling the present time (z — 0) as,A value to that rele- 
vant to each DE model by using eq.(6.3) at the present 
time. We again have that: 



08, DE = 08, A' 



^DE.O) P DE ,0 *(On,i?8 



Da(0) I/ ^A,0 *(-2 m ,A,^8 



(6.8) 



with f2 TO ,A denoting the value for the reference A model 
(in our case fl m ,A — 0.28, see section 3), while we have 
(for fixed Hq, as in our case), that: 



-PpE.o 



1/2 



n 



m,A 



-2, 



(6.9) 



It thus follows that if the o^a value is known, one may de- 
rive the corresponding <78,de values to the other DE mod- 
els. In this concern, the combined SNIa+BAO+WMAP5 
analysis of Komatsu et al. [5] (see also [94]) pro- 
vided a value of a$,A — 0.812, while the correspond- 
ing WMAP7 analysis yielded [for w(z) = —1]: er 8 ,A — 
0.803 (using the WMAP7 alone) and a s ^ 0.807 for 
the joint WMAP7+BAO+i? analysis [5]. A recent 
analysis based on cluster abundances have also fur- 
nished the following degenerate combination: <78,a = 
0.83 ± 0.03(r2 m , A /0.25)-°- 41 [95], which for our 'case 
(ft mA = 0.28) implies cr 8 ,A ^ 0.792. Fu et al. [96] 
based on a weak-lensing procedure found cts.a = 0.837 ± 
0.084(ft m:A /0.25)~ a53 which implies <t 8 ,a ~ 0.788 in our 
case. In addition, studies based on the peculiar ve- 
locities statistical analysis [97] obtained <T8,a = 0.88 ± 



0.05(-2 TO) a/0.25) 



-0.53 



or (is, a — 0.829 in our case. It 



should be stressed that the the average scattering of the 
four independent as, a values (of those based on WMAP, 
we use only the joint WMAP7+BAO+ff result) is quite 
small (((78. a) — 0.804 ± 0.018), thereby reinforcing the 
consistency of the different measurements. Finally, by 
inserting the latter value in Eq.(6.8) we can estimate the 
corresponding cts.de values, listed in Table 1, to be used 
in our mass function analysis. 

For completeness, it should be remarked that some re- 
cent analyzes are suggesting significant higher values of 
ag- For example, Watkins et al. [98] studying the bulk 
flow on scales of <~ 100 /i _1 Mpc found a eg normaliza- 
tion which is increased by a factor of <~ 2 with respect 
to the one of ACDM model. If these results are correct, 
the ACDM model would be strongly challenged. How- 
ever, such a discussion is beyond the scope of the current 
work. 



B. Halo Mass Function & Number Counts of DE 
Models 

We now pass to our results. In the upper panels of 
Fig. 1, for two different rcdshifts (z = 0.5 and 2), we 
display the integral halo mass function, n(> M), for all 
DE models previously discussed. The different models 
are characterized by the symbols and line- types presented 
in Table 1. 

In the corresponding lower panels one may see the frac- 
tional difference between each DE model with the concor- 
dance ACDM model, that is, 5n/n\ — (tide — ^a)/«a- 
We stress that we have shown the case z — 2 only for 
comparison of the model expectations. In other words, 
it is not a statement of the viability to actually observe 
clusters at such a large redshift (not only due to tech- 
nical limitations but also because the majority of dark 
matter halos are not expected to have virialized at such 
a redshift, thereby producing their X-ray signatures). 

It is worth noticing the mass function expectation of 
the Aps model (open red circles), the BRG model (blue 
dashed line) and, at high redshifts, of the CCDM model 
(magenta stars) appear to be completely different from 
the reference ACDM cosmology. 

Given the halo mass function we can now derive an 
observable quantity which is the redshift distribution of 
clusters, Af(z), within some determined mass range, say 
Mi < M/Mq < M 2 . This can be estimated by integrat- 
ing, in mass, the expected differential halo mass function, 
n(M, z), according to: 



dV f M2 

Af(z) = — / n(M,z)dM, (6.10) 
dz J Ml 

where dV/dz is the comoving volume element, which in 
a flat universe takes the form: 



with r(z) denoting the comoving radial distance out to 
redshift z: 



dV . 2/ \dr 
-=4nr (z)-(z), 



r(z) 



c 

Hi 



o Jo 



dx 



(6.12) 



In Fig. 2 (upper panel), we show the theoretically ex- 
pected cluster redshift distribution, Af(z), for all the 
models studied (line and symbol types are listed in Table 
1), for cluster of galaxies size halos, ie., Mi = 10 14 M Q 
and M 2 = 10 16 M . It is evident that many of the 
different models show significant differences with respect 
to the concordance A model and therefore we could, in 
principle, distinguish them. In the lower panel of Fig. 2 
we show the relative differences of the various DE mod- 
els with respect to the expectations of the concordance 
ACDM model. Note that the three models (A PS , CCDM 
and BRG models), for which we found very large mass 
function differences (Fig. 1), are not shown because their 
large relative differences are beyond the limits of the 
lower panel. 

Let us now discuss the expected redshift distributions, 
based on two future cluster surveys, and also the pos- 
sibility to observationally discriminate the different DE 
models. 

These two realistic future surveys are: 

(a) the eRDSITA satellite X-ray survey, with a flux limit 
of: /iim = 3.3 x 10~ 14 ergs s _1 cm~ 2 , at the energy band 
0.5-5 keV and covering ~ 20000 deg 2 of the sky, 

(b) the South Pole Telescope SZ survey, with a limiting 
flux density at vq = 150 GHz of / yo ji m = 5 mJy and a 
sky coverage of ~ 4000 deg 2 . 

To realize the predictions of the first survey we use the 
relation between halo mass and bolometric X-ray lumi- 
nosity, as a function of redshift, provided in [100], ie: 



L{M, z) = 3.087 x 10 4 



ME(z) 



1.554 



ergs 



lO^/i-iMg 

(6.13) 

The limiting halo mass that can be observed at redshift z 
is then found by inserting in the above equation the lim- 
iting luminosity, given by: L = 47rd 2 fn m Cb, with dh the 
luminosity distance corresponding to the redshift z and q, 
the band correction, necessary to convert the bolometric 
luminosity of eq.(6.13) to the 0.5-5 keV band of eROSITA. 
We estimate this correction by assuming a Raymond- 
Smith (1977) plasma model with a metallicity of O.4Z , 
a typical cluster temperature of <~ 4 keV and a Galactic 
absorption column density of uh = 10 21 cm~ 2 . 

The predictions of the second survey can be realized 
using again the relation between limiting flux and halo 
mass from [100]: 



/i/o,lim 



2.592 x 10 8 mJy 



M 



E 2 l z (z) 



(6.11) 



d\{ z ) Vio 15 m , 

(6.14) 

where dA(z) = g?l/(1 + z) 2 is the angular diameter dis- 
tance out to redshift z. 









1000 - 

100 

10 
0.4 

< 0.2 

> 



-' ' 1 1 1 1 1 ' " 1 1 1 1 1 1 1 1 

: / z 

-_ y z 




: * * „ x : 

-, , , , , I , , , , I , , ,- 


\ u □ * * x J - 



0.5 1 1.5 0.5 1 1.5 
z z 

FIG. 3: The expected cluster redshift distribution of the dif- 
ferent DE models and for the two different future cluster sur- 
veys (upper panels), and the corresponding fractional differ- 
ence with respect to the reference ACDM model (lower pan- 
els). Symbols correspond to the different DE models as indi- 
cated in Table 1. 



differences of the DE models but only up to a redshift at 
which they are significant, that is, such that: 



A/de - A/a 



(A/de -A/a) 1 / 1 



> 3.5, 



(6.15) 



FIG. 2: The expected redshift distribution of M<; 10 14 M Q 
clusters (upper panel) of the different DE models and the cor- 
responding fractional difference between the models and the 
reference ACDM model (lower panel). The lower panel shows 
only those DE models that have fractional relative differences, 
with respect to the ACDM model, of <j 45%. Symbols corre- 
spond to the different DE models as indicated in Table 1. 



In Fig. 3 (upper panel) we present the expected red- 
shift distributions above a limiting halo mass, which is 
Mi = A/ii m it = max[lO 14 M0, Mf], with Mf correspond- 
ing to the mass related to the flux-limit at the different 
redshifts, estimated by solving eq.(6.13) and eq.(6.14) for 
M. In the lower panels we present the fractional differ- 
ence between the different DE models and ACDM, simi- 
larly to Fig. 2, but now for the realistic case of the previ- 
ously mentioned future cluster surveys. It is evident that 
the imposed flux-limits together with the scarcity of high- 
mass halos at large redshifts, induces an abrupt decline 
of Af(z) with z, especially in the case of the eROSITA X- 
ray survey (note the shallower redshifts depicted in Fig. 3 
with respect to Fig. 2). 

In the lower panels of Fig. 3 we display the relative 



with Ade the redshift distribution predicted by some 
DE model and A/a the corresponding redshift distribu- 
tion of ACDM model. However, such a criterion does 
not take into account cosmic variance and possible ob- 
servational systematic uncertainties which can hamper 
detecting small (but according to Eq. 6.15 significant) 
relative differences. We believe that relative differences 
of <S 5% will be difficult to detect especially at relative 
high redshifts. 

In Tabic 1, one may see a more compact presentation 
of our results including the relative fractional difference 
between all DE models and the ACDM model, in 3 dis- 
tinct redshift bins and for both future surveys. 

Based on our Af(z) analysis and the results presented 
in Figure 3 and Table 1, we can now divide the studied 
DE models into those that can be distinguished obser- 
vationally and those that are practically indistinguish- 
able from the current paradigm (ACDM). The latter DE 
models are the following four: XCDM, HRDE, TCM 
and PNGB. One has to remember however that these 
results are based on using DE model parameters that 
have been fitted by the present day cosmological data 
(see section 4). As an illustrative example, we remind 
the reader that the XCDM model, compared here, is one 



with w = —0.99; if future cosmological data would pro- 
vide a different value for the equation of state parameter 
then the N(z) predictions of such an XCDM model could 
be quite different than those derived here. 

Regarding the models that are distinguishable with re- 
spect to the concordance model, three of them (BRG, 
Aps and CCDM) show extremely large variations mak- 
ing it trivial to distinguish them. From the rest of the 
distinguishable DE models all of them show clear signs of 
difference, at all redshifts and in both future cluster sur- 
veys, with respect to the ACDM expectations; only the 
Arg model needs to be distinguished at higher redshifts 
0.5). 

As an additional test and in order to check the sensi- 
tivity of our results only on the functional form of the DE 
equation of state parameter, we have imposed a unique 
value of the rms mass fluctuation normalization to all 
DE models according to: (j&,de = &8,A, and repeated 
the cluster-size halo N(z) analysis. We now find slightly 
different results although in the same overall direction 
with our main analysis. For example, the DE models 
that cannot be distinguished from the reference ACDM 
model are now the XCDM, TCM, PNGB, EDE and A RG , 
with the first three being common in both analyzes, the 
Afic model showing a small difference and the EDE be- 
ing the only model that shows a significantly different 
behavior between the two analyzes. 

Finally we would like to mention that an interesting 
paper appeared recently [101] and among other issues, it 
compares different forms of the halo mass function and 
its redshift evolution using N-body simulations of the 
ACDM and wCBM (to =const) models. They do find 
some differences at the few percent level. Although our 
analysis is self-consistent, in the sense that we compare 
the expectations of DE models with respect to those of 
the concordance cosmology using the same mass function 
model, we plan to investigate in a forthcoming paper how 
sensitive are our observational predictions to the different 
mass functions fitting formulas. 



obtaining tight constraints on the main cosmological and 
dark energy free parameters. At the level of the resulting 
Hubble function, we have found that of all dark energy 
models (apart a Brane world cosmology), are statistically 
indistinguishable (within la) from a flat ACDM model, 
as long as they are confronted with the quoted set of 
observations. 

On the other hand, despite the fact that these models 
closely reproduce the ACDM Hubble expansion, we show 
that eight out of the twelve studied DE models, using the 
observationally fitted model parameters, can be differen- 
tiated from the reference ACDM model on the basis of 
their redshift distribution of cluster-size halos. Such a 
comparison was made possible by using the expectations 
of a future X-ray (based on the eROSITA Satellite) and 
SZ cluster surveys (based on the South Pole Telescope). 

The main comparison results can be summarized in the 
following statements (for nomenclature see section 4): 

• Four DE models, namely, XCDM, HRDE, TCM 
and PNGB cannot be distinguished from the 
ACDM model at any significant level. 

• Seven models, ie., the BRG, CPL, LRDE, EDE, 
VCG, Aps, CCDM, can be easily distinguished due 
to the fact that they show strong and significant 
variations with respect to the concordance A model 
even at z = 0, implying that even with the present 
day surveys one could effectively distinguish them. 

• The Ajig model, although presents relatively small 
variations with respect to the concordance A model, 
it can be clearly distinguished at relatively high 
redshifts (z^ 0.5). 

In a future work we will present a comparison be- 
tween our model predictions and the observationally de- 
termined cluster mass function at different redshifts as 
well as the available (X-ray or optical) cluster redshift 
distribution. 



7. CONCLUDING REMARKS 

In this paper, we have investigated the cluster abun- 
dances beyond the conventional ACDM cosmology by us- 
ing several parameterizations for the dark energy. In 
order to do that, we first performed a joint likelihood 
analysis using the most recent high quality cosmological 
data (SNIa, CMB shift parameter and BAOs), thereby 
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TABLE I: Numerical results. The I s * column indicates the DE model. 2 nd column are symbols or line types of the models 
appearing in Figs. 2, 3 and 4. 3 rd and 4 th columns show the os.de and <5 C values. The remaining columns present the fractional 
relative difference between the models and the ACDM cosmology for two future cluster surveys discussed in the text. Bold 
letters denote those models that can be clearly distinguished from flat ACDM model at some redshift bin (uncertainties appear 
only when their value is > 10~ 2 ). 



